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Counting dimensions 
of L-harmonic functions 

O ■ By Peter Li and Jiaping Wang* 
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In this article, we will consider second order uniformly elliptic operators 
of divergence form defined on R" with measurable coefficients. Mainly, we 
will give estimates on the dimension of space of solutions that grow at most 
polynomially of degree d. More precisely, in terms of a rectangular coordinate 
system {xi, . . . ,Xn}, a second order uniformly elliptic operator of divergence 
form, L, acting on a function / E H^^^{R^) is given by 



n 



CN ' where (a*-'(x)) is an (n x n)-symmetric matrix satisfying the ellipticity bounds 

O. (0.1) A/<(a*^)<A/ 

^^ I for some constants < A < A < cxd. Other than the ellipticity bounds, we only 

T^ ' assume that the coefficients (ay) are merely measurable functions. 

C^ ■ A function / G H^^^(R") is said to be L-harmonic if Lf = 0. Since (aij) 

^ , are only measurable, this equation is being interpreted in the weak sense. The 

celebrated De Giorgi-Moser-Nash theory asserts that any weak solution / must 

be C" for some < a < 1. A global version of this theory implies that there 

exists < a < 1 such that if p{x) denotes the Euclidean distance from x to 

the origin, then any L-harmonic function / satisfying the growth condition 

|/(x)|=0(p"(x)). 
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as X ^ oo, must be identically constant. This statement can be thought of as 
a C" regularity result at infinity. In general, for each real number d > 0, we 
can define 

Tid^L) = {/ G i^iL(K") I ^/ = and |/(x)| = 0(p^(x))} 

to be the space of all polynomial growth L-harmonic functions of degree at 
most d. Let us denote the dimension of TidiL) by hd{L). The above growth 
result can be restated as 

hd{L) = 1 

for all < d < a < 1. The main purpose of this article is to give an upper 
bound estimate on hd{L). 

Before stating our result precisely, let us give a brief outline on the his- 
tory of this problem. After more than 20 years since the establishment of the 
De Giorgi-Moser-Nash theory, Avellaneda and Lin [A-Ln] studied the solution 
space of a uniformly elliptic operator of divergence form. They assumed that 
the coefficients are Lipschitz continuous and periodic. Using the theory of ho- 
mogenization, they showed that every polynomial growth L-harmonic function 
is necessarily a polynomial of periodic coefficients. Moreover, there is a linear 
isomorphism between the harmonic polynomials and the polynomial growth 
L-harmonic functions on R". In particular, this implies that hd{L) are finite 
for all d. Later, Moser and Struwe [M-S] using a different argument proved the 
same result without the Lipschitz continuity assumption on the coefficients (see 
also [Ln] for this last point). Their argument also works for a certain class of 
nonlinear operators. Recently, based on the concept of G-convergence, Lin [Ln] 
considered the case when L is an asymptotically conic operator. In this case, 
his results imply that hd{L) are finite and have an upper bound depending on 
n, d, A and A. Zhang [Z] followed this direction by assuming that L converges 
to a possibly S-*^ family of conic operators and established some similar results. 
Last year, the first author [LI] proved that for a general uniformly elliptic 
operator of divergence (and nondivergence form) with measurable coefficients, 
there exists a constant C > depending only on n, A and A, such that, 

(0.2) hd{L)<Cd''-^ 

for ah d>l. 

Note that when L is the Euclidean Laplacian on R", one can compute 
hd{A) directly because 7id{A) is spanned by the homogeneous harmonic poly- 
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nomials of degree at most d. In fact, it is known that 
fn+d-l\ (n+d-2 

= ^^^^-^id+n-2)id + n-3)...id+l) 
[n 1). 

^ (d"-i + 0{d''~^)) . 



{n-l)\ 



This example indicates that the estimate (0.2) is sharp in terms of the order 
of d. This computation also raises the question of whether it is possible to 
obtain an estimate for h(i{L) which will reflect the sharp constant. One may 
ask if there exists a constant C{n, A A~^) > depending only on the prescribed 
quantities and has the property that C{n, A A~^) ^ 1 as A A~^ -^ 1, such that, 

/irf(L) < C(n,AA-i) ^^— d"-\ 
[n — 1)\ 

In particular, such an estimate will be sharp as realized by L = A. 

This article is an attempt to prove the validity of this estimate. While we 
can not establish the estimate in the precise form, we do prove a sharp upper 
bound for the truncated sum of hd{L). In fact, our estimate only depends on 
the ratio of the ellipticity constants at infinity. More precisely, let us denote 
Ar and A^. to be the ellipticity bounds satisfying 

\rl< {a'\x)) <ArI 

for all X G W^ \ B{r) with B{r) = {x \ p{x) < r} being the Euclidean ball 
of radius r. Obviously, A,, and A^ are montonicially increasing and decreasing 
functions of r, respectively. Also, they satisfy 

< A < A,. < Ar < A < cx). 

Let us define their limits at infinity by 

Aoo = lim Ar, 

1 — >oo 

and 

Aoo = lim ^r- 

r— »oo 

The main result of this article is to prove that 

i— 1 

This upper bound is sharp for L = A by simply summing up the formula for 
hdiA). 
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We would like to point out that a similar type of estimate was proved by 
the authors [L-W2] for the Laplacian on a complete manifold with nonnegative 
sectional curvature. In both the nonnegatively curved manifold case and our 
present case, our methods rely on restricting the operator on codimension-1 
submanifolds. A major difficulty for measurable coefficients is that the coeffi- 
cients might not be measurable when restricted on a codimension-1 submani- 
fold. As a result, extra work is needed to overcome this problem. In a parallel 
manner, the lack of regularity issue occurs on the cut-locus for the manifold 
case and extra but different effort was required to overcome that hurdle. 

It is also worth mentioning that tremendous progress has been made in 
the direction of the establishment of Liouville properties and in counting the 
dimensions of polynomial growth harmonic functions on a manifold. We will 
refer the interested reader to the paper [C-M] by Colding and Minicozzi and 
a recent survey article by the first author [L2], and some subsequent articles 
[S-T-W], [L-Wl]. 

1. (/>-harmonic functions on a complete manifold 

Our strategy is to translate the problem of studying L-harmonic func- 
tions on M" to that of studying i;^-harmonic functions on a manifold (M",(7). 
Let (M",5) be an n-dimensional Riemannian manifold. Let (/> be a positive 
measurable function on M. A function / is said to be (/)-harmonic if it satisfies 

div(0V/) = O, 

where V and div are the gradient and the divergence with respect to the Rie- 
mannian metric g. In terms of local coordinates {xi, . . . , a;„}, the Riemannian 
metric can be expressed as g = Y^^ j=i 9ij dxi dxj. Using this, we can write 

(1.1) av„v/)=t;^|-(*v^»-^||). 

where G = det (gij) and {g"-^) = {gij)~^ . 

For any nonnegative real number d, we denote 

■HdigA) = {/|div(</.V/) = 0, |/(x)| = 0{r\x))} 

to be the space of (;^-harmonic functions with at most polynomial growth of 
degree d with respect to the geodesic distance r(x) to some fixed point. Let 
us denote the dimension of TCdig, 4') by hd{g, (p) = dimTidig, (j))- 

To set up the L-harmonic problem as a 0-harmonic problem, we recall 
that p{x) is the Euclidean distance function to the origin. Let us introduce the 
Riemannian metric g = Y17j=i 9ij ^-^i ^-^j ^^ ^'^ by 

(1-2) {g^^) = {g,,r'=w{x){a^n, 
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where 




(1.3) 


, , _ ■^ ij dp dp 
^ dxi dxi' 


If we define 




(1.4) 


<f>g{x) = w 2 (x) Jdet{a'J), 



then one easily verifies by (1.1) that / is L-harmonic if and only if it is 
(/ig-harmonic on M = (M",^). Since the metric g is uniformly equivalent to 
the Euclidean metric, we have 7id{L) = Ti-dig, 4'g) for all d > 0. Our aim is to 
work with hd{g,(j)g) instead. Let us first establish the following lemma. 

Lemma 1. Let T a Q he a smooth subdomain properly contained in an 
open domain fi C M. Let (j) he a positive measurable function whose restriction 
to the boundary dV of T is also measurable. Suppose H is a k- dimensional 
space of nonconstant cp-harmonic functions defined on Q, and 

= r]i <r]2 <r]3 < ... 

are the eigenvalues of dT with respect to the operator (/>~^ div (</> V), where V 
is the tangential gradient on dV. Then for any orthonormal basis {ui, . . . , u^} 
of H with respect to the inner product 

D{u,v) = I {Vu,Vv)(j)dV, 

we have 

k 



rr T^ JdV 



where dA is the area element on dV induced by the metric g. 

Proof. Let us first point out that the right-hand side of the inequality is in- 
dependent of the choice of the D-orthonormal basis. Let {wi, . . . ,Wk^i} be the 
first k — 1 eigenfunctions on dV corresponding to the eigenvalues {771, . . . , rjk-i}. 
Using a similar argument as in Lemma 1.2 of [L-W2], we can find a 
D-orthonormal basis {ni, . . . , u^} of iif such that 

Ui Wj 4> dA = 
lav 

for all 1 < J < i < A;. Hence by the variational principle, we obtain that 

(1.5) r]i / u^i(t)dA< / \Vui\'^(t)dA. 

Jar Jar 
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On the other hand, we have 

k k 

k 

Applying (1.5) to the first term on the right-hand side and then combining the 
two terms, we conclude that 

2V??f<V/ \Vui\'^(l)dA. D 



2. L-harmonic functions 

We are now ready to prove the main theorem. The manifold M = (R", g) 
is assumed to be the R" endowed with the Riemannian metric satisfying (1.2) 
and (1.3). The function (pg is defined by (1.4). 

Theorem 2. Let h^ = dmi7id{g,4'g) be the dimension of the space of 
polynomial growth (/)g-harmonic functions of at most degree d. Then for any 
sequence of nonnegative numbers = gq < ai < • • • < aj = d, ha must satisfy 

5^(a. - a._i) ha^_, < h^ - (d + 2n - 1)«. 

In particular^ when d is a positive integer, by taking Oj = i, we see that 






Proof. Let us define 

Ki9,4>9) = {f^nd{g,<t>g)\f{o) = o} 

and its dimension is denoted by h'^ = diniH'^{g,(l)g). Obviously h'^ = h^ — 1, 
and we will estimate h'^ instead. Let us denote the Dirichlet form on B{r) by 



Dr{flj2)= [ {Vgfl,Vgf2)^gdVg. 

JB{r) 
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Using the sequence = ao < ai < 02 • • • < flj = d, we can decompose the 
space ^{'^{g, cpg) with respect to the inner product Di into a direct sum 

where each 

Hi = Ha^ig, 4>g) Q na,_i{g, 4>g) 

is the subspace consisting of harmonic functions of growth order between Oj-i 
and Oj. Let us denote its dimension by dimi^j = ki. If we view the Dirichlet 
form Dr as a bihnear form defined on the space 'H'^{g, (j)g), then the determinant 
of Dr when computed with respect to a Di-orthonormal basis must satisfy the 
growth assumption 

(2.1) detDi^r <Cr'', 

where s = X]i=i(2(aj - 1) + n) ki. 

For a fixed r > 0, and for any e > 0, there exists a smooth matrix {¥^) 
and a set A C B(r) with its measure satisfying \A\ < e such that 

A/< (&*■'■) < A/ 

on B{r), and 

\a'^ -b'^\ <e 

on B{r) \ A. Moreover, for any fixed vq < r, we may also insist that 

onB{r)\B{ro). 

Let us define the corresponding Riemannian metric h = Y17i=i ^*i '^•^« ^-^j 
on R" by 



with 



n 



Wh{x) = Y, b'' 



dp dp 

dxi dxi 



We also define the corresponding function 



n-2 I 

(l)h{x)=Wf^^ (x)^det(6*J). 

Note that we have 

A < ty/i < A 

and 
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For each u G 7{'^{g,(f)g), let v G H^{B{r)) be the solution to the Dirichlet 
problem 



(2.2) 



i^(-<o- 



on B{r) and 

(2.3) v = u 

on dB{r). We denote /C(r) to be the space consisting of all such solutions v. 
The maximum principle then asserts that dim/C(r) = h'^. For tq < t < r, 
applying Lemma 1 to H = ]C{r) and F = i?(t), we get 



d -, dp 

(2.4) 2Vr??(e,t)< V / |V;,n,|2 </.;, d^, 

where ryj(e, t) are the eigenvalues of the operator cftf^ div((?!>fe V^) on dB{t) with 
V/i being the tangential gradient of the metric h on dB{t). 

If we define the Dirichlet form on the space /C(r) with respect to h by 

Dt{h){fij2)= [ {Vhfi,Vhf2)cl>hdVh, 

JB{t) 
then 

K 

4>hdAh. 



a /» 

(lndet^,(,)A(/i))' = j;/ 



r[JdB{t) \^hp\ 

Combining with (2.4) and the fact that 

— ^Vo' 
we conclude that 

2^r/f (e,t) < A,„ (lndeto^(;,) A(/j))'. 
Integrating from rg to r, we get 

K 

(2.5) 2 rj^r]fie,t)dt < A,„ lndet^^^(^)D,(/i). 

We now claim that as e ^ 0, the solution v G H^{B{r)) to the boundary 
value problem (2.2) and (2.3) must satisfy 

(2.6) / \Vg{u-v)\'^dVg^O. 

JB(r) 
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Indeed, 



n It- 

/ y v= 



jd{u — v) d{u — v) 



dxi dxi 



dVo 



-X„<— '1.1- (<--"' J^)-' 



where dVo is the volume form with respect to the standard Euchdean metric. 
Sphtting this integral over the two domains B[r) \A and A, we see that there 
exists a constant 5 > with (5 ^ as e ^ 0, such that, 






<e / \Vo{u-v)\\VQu\dVo 

JB(r) 



1 1 

2 / n \ 2 

|2 JT/ 1 I / IV7 „,|2 



and 



<e { \Vo{u-vrdVo\ / \Vou\'dVo\ , 

lB{r) J \JB{r) I 



<C \Vo{u-v)\\Vou\dVo 

J A 

<clj \V^{u-v)\^dvA (f\Vou\^dVo^' 



<C6 I \Voiu-v)\'dVo 

\JB(r) . 
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Using the fact that the metrics g and h are uniformly equivalent to the Eu- 
clidean metric, we conclude that 



\Vg{u-v)\UVg<C \Vh{u-v)\UVh 

B{r) JB{r) 



1 

\ 2 
2 



<C6 { \Vo{u-v)\^dVo 

'B(r) 



1 
\ 2 
2 



<C5\ \Vg{U-V)\'dVg\ , 

\JB{r) J 

and the claim is proved. In particular, we conclude from (2.6) that 
(2.7) / iVhvfdVh^ f \VgUfdVg 



JB(t) JB{t) 

as e ^ 0. 

We now claim that for each A; > 1, the k^^ eigenvalue rik{e,t) must satisfy 

(2.8) rik{e,t)>Xl^^r,k{l)t'\ 

where ryfc(l) denotes the k Euclidean eigenvalue of dB{l) with respect to the 
induced Euclidean metric. In fact, the variational characterization implies 

rik[e,t) = mt sup 



w f(zw infceR /g^(j) (/ - c)2 (j)h dAh '' 



where the infimum is taken over all /c-dimensional subspace W in the space of 
smooth functions on dB{t). By the co-area formula, we have 



\Vhp\'^ dAh = y det(/iij) dAo, 

where dA^ and dA^ are the area elements on the sphere induced by the metric 
h and the Euclidean metric respectively. Thus, by the choice of (ph) we have 



So, 



r]k{e,t) = inf sup 



dAo = 4>hdAh. 

IdB(t)\^hf\'dAo 



w f^w infceR fdBit)(f ~ c)^ "^^0 



''" V f(,w infceR fdB{t) (/ ~ ^)^ ^^0 



>^ro Vk{t) 



A' %(l)t-^ 
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where ??fc(t) is the k EucHdean eigenvalue of the sphere dB{t). This proves 
the assertion (2.8). 

In is known that [B-G-M]) for each A; > 1, 

Vk{l) = q^ + {n - 2)q 

where q is the least positive integer satisfying 

However, one checks easily that 

2 / n-1^"-^ 



hence we have 



and 



Therefore, 



1 



'in- 1)! ,\ "-1 n- 1 



,1(1). (^^•) 



2 ' 
(n — 1)! \ "^ n — 1 



(2.10) ^,4a)>((!i^)-''^,:i-!i^t 

,(n-l)!\^ /■'^-i ^ , n-1 , 

- 2 ; io 2 

1 
(n— l)!\"-i n— 1,, ,^i_ n— 1, 
' (/c - l)"-i fc 



2 y n ' ' 2 

, (n — 1)!\ "-1 n — 1 , ^j_ , , , 

>(^^j __/c.-i-(n-l)fc. 



Applying (2.8), (2.10) to (2.5), we have 



(2.11) 2A,, ni!i^y" !L_l^h',)^^.^n-l)h'A (Inr-lnro) 






<2/ j;r/?(e,t)(it 
< A^y lndet£,^^^(;,)Z)^(/i). 
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Letting e ^ and using (2.7), we conclude from (2.11) that 

< (j2{2{a^-l)+n)kA (Inr + C), 

where we have used (2.1) and the constant C is independent of r. Dividing 
both sides by Inr and letting first r ^ oo and then tq ^ oo, we get 



["[ . "] ' n)ki. 
1=1 



< 5^(2(0,-1)+^ 



Using the fact that fcj = h'^. — h'a^_^, we can rewrite this inequality into the 
form 



(2.12) Y.^a, - a,_i) h',^_^ < U - l + \ h'^ 



i=l 

1 



(^) (^)^^<^^)^-<"-"< 



^(-T-)^^-fe)(^)^^<^^ 



d 

1 

i" n - 1 ,. , 

n-l 



As a function of /i^, the right-hand side of (2.12) achieves its maximum when 



h'^ satisfies 



^oo 



1 



^-^ ^^"-'^'V-(;,.);^_r^+^_2)=0. 



A^/ V 2 / ^ "' V 2 



Plugging this value of h', into (2.12) and simplifying, we obtain 



|:«.,-<..-)'c.4(fe)""(''-f-^)' 
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Note that h'^ = hd — 1; thus we have 
j 3 



i=l i=l 

n— 1 



This proved the first assertion. In the case d is a positive integer and Oj = i, 
the first assertion implies that 

n— 1 

2n - ir. 



1=1 



Now replacing d by d + 1, we get the second assertion. D 

Corollary 3. Let 






be a uniformly elliptic operator defined on R" with measurable coefficients sat- 
isfying (0.1). // hd{L) = di\ra.'Hd{L) denotes the dimension of the space of 
polynomial growth L-harmonic functions of at most degree d. Then 

d / A \ n-l 



j=l 
and 



T.u^^^{tT ^(^"+^(^"")) 



hminf d-("-i) /i,(L) < i^X ' -^ 
d^oo "'^ ^ - VAoo/ (n-1)! 
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